arXiv: 1503.00513v2 [hep-ph] 29 May 2015 


Factorization model for distributions of quarks in hadrons 

B.I. Ermolaev 

Ioffe Physico-Technical Institute, 194021 St.Petersburg, Russia 
M. Greco 

Department of Mathematics and Physics and INFN, University Roma Tre, Rome, Italy 

S.I. Troyan 

St. Petersburg Institute of Nuclear Physics, 188300 Gatchina, Russia 

We consider distributions of unpolarized (polarized) quarks in unpolarized (polarized) hadrons. 
Our approach is based on QCD factorization. We begin with study of Basic factorization for 
the parton-hadron scattering amplitudes in the forward kinematics and suggest a model for non- 
perturbative contributions to such amplitudes. This model is based on the simple observation: after 
emitting an active quark by the initial hadron, the remaining set of quarks and gluons becomes 
unstable, so description of this colored state can approximately be done in terms of resonances, 
which leads to expressions of the Breit-Wigner type. Then we reduce these formulae to obtain 
explicit expressions for the quark-hadron scattering amplitudes and quark distributions in Kt- and 
Collinear factorizations. 

PACS numbers: 12.38.Cy 


2 


I. INTRODUCTION 

QCD factorization, i.e. separation of perturbative and non-perturbative QCD contributions, proved to be an efficient 
instrument for describing hadron reaction at high energies. Being first applied to processes in the hard kinematics in 
the form of Collinear factorizationjlj, it was soon extended to cover the forward kinematic region, with DGLAPfl] 
used to account for perturbative contributions. Then, in order to be able to use BFKLQ, a new kind of factorization, 
Kt -factorization was suggested in Ref.Q. These kinds of factorization are usually illustrated by identical pictures. 
For instance, factorization of the DIS hadronic tensor W^ is conventionally depicted by the construction in Fig. |T] 
both in Collinear and in Kt- factorizations, where the upper, perturbative blob and the lower, non-perturbative blob 
are connected by two-parton state. The upper blob in Fig. Q] is calculated with regular perturbative means. On the 





FIG. 1. Conventional illustration of QCD factorization. The s-cut of the graph is implied. 

contrary, the lower blob is conventionally introduced from purely phenomenological considerations. Collinear and 
Kt- factorizations operate with different parametrizations for momentum k of the connecting partons and as a result, 
they are described by different formulae. Collinear factorization assumes that 


k = Pp, 

while Kt -factorization allows for the transverse momentum in addition: 


(1) 


k = /3p+ k±, (2) 

accounting therefore for one longitudinal and two transverse components of k. However as a matter-of-fact, k has four 
components: two of them are longitudinal and the other two are transverse. Accounting for the missing longitudinal 
component a (for definition of a see Eq. (|3|l) drove us to suggesting a new, more general factorization which we 
named in Ref. [|| Basic factorization. In contrast to Kt~ and Collinear factorizations, the analytic expressions in 
Basic factorization can be obtained from the graphs of the type of the one in Fig. [T| with applying the standard 
Feynman rules. 

It is worth reminding briefly our derivation of Basic factorization, for detail see Ref. Q. Let us consider the 
Compton scattering amplitude off a hadron in the forward kinematics. It is depicted in Fig. [2] 



P P 


FIG. 2. Amplitude for forward Compton scattering off a hadron target. 

The blob in Fig. [2] denotes ensemble of perturbative and non-perturbative contributions. This blob can be expanded 
into an infinite series of terms, each of them is represented by two blobs connected with n parton lines, n = 2,3 ,... 
Considering only the simplest, two-parton state, we arrive to the graph similar to the one in r.h.s of Fig.|T]but without 
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the s-cut and with the both blobs accommodating perturbative and non-perturbative contributions at the same time. 
The integration of the convolution in Fig. [T] over momentum k now runs over the whole phase space and it is expected 
to bring a finite result. However, the propagators of the connecting partons become singular at k 2 = 0 (we neglect 
quark masses). Besides, the upper blob may contain IR-sensitive perturbative contributions ~ ln™(2 pk/k 2 ) (with 
n = 1,2,..). In addition, it yields the factor 2 qk/k 2 , when unpolarized gluon ladders are included into consideration. 
The only way to kill such IR singularity is to assume that the lowest, non-perturbative blob should tend to zero fast 
enough when k 2 —> 0 . Doing so and repeating a similar procedure to regulate the UV singularity, we bring the 
convolution in Fig. |T] to agreement with the factorization concept: perturbative and non-perturbative contributions 
are located indifferent blobs. This is a new form of QCD factorization which we name Basic factorization. 

We demonstrated in Ref. [EJ that Basic factorization can be reduced step-by-step first to Kt- and then to Collinear 
factorizations. In Ref. (5| we began with considering Basic factorization for Compton scattering amplitudes in the 
forward kinematics, where integration over momentum k of the connecting partons in Fig.[l]runs over the whole phase 
space. Confronting two obvious facts that, on one hand, the integration over k should yield a finite results and that, 
on the other hand, the perturbative part in Fig. |T| (the upper, perturbative blob and propagators of the connecting 
partons) is divergent in both the infra-red (IR) and ultra-violet (UV) regions, allowed us to impose integrability 
restrictions on the lowest blob, which are necessary for the convolution in Fig. |T] to be finite. The obtained restrictions 
led us to theoretical constraints on the fits for the parton distributions to the DIS structure functions in Collinear 
and Kt~ factorizations. In particular, we predicted the general form of the fits in Kt -factorization and excluded the 
factors x~ a from the fits in both Kt- and Collinear factorizations. 

Another interesting object, where factorization is used, is distributions of partons in hadrons. In the present paper 
we examine their properties in IR and UV regions and suggest a simple resonance model for the non-perturbative 
contributions to the parton distributions. Our argumentation in favor of this model is as follows: after emitting an 
active quark by a hadron, the remains of the hadron, i.e. a set of quarks and gluons, acquires a color and therefore 
it becomes unstable. So, this colored state can be described in terms of resonances. We begin with considering 
amplitudes of the quark-hadron (QHA) and gluon-hadron (GHA) scattering in the forward kinematics. The Optical 
theorem relates such amplitudes to the parton distributions. Throughout the paper we use the standard Sudakov 
parametrizationQ for momentum k of the connecting partons: 

k = —aq' + flp' + k_ l, (3) 

where momenta q' and p' are massless, p' 2 ~ q' 2 ~ 0, and they are made of the hadron momentum p and the parton 
momentum q: 

p' =P + q, q' = q + xip, (4) 

where x 2 = —p 2 /w = — M 2 /w , x\ = —q 2 /w, with w = 2 pq « 2 p'q'. In these terms 

2pk = w{—a — £2/3), 2qk = w(/3 — x±a), k 2 = —wa/3 — k\. (5) 

In Sect. II we introduce the quark-hadron scattering amplitudes in the forward kinematics and examine their IR and 
UV behavior. In Sect. Ill we consider separately the unpolarized and spin-dependent quark-hadron amplitudes in Basic 
factorization and suggest a model for non-pertubative contributions to the amplitudes. This model involves a spinor 
structure accompanied by invariant amplitudes and T^ s \ In Sect. Ill we specify the spinor structure of the non- 
perturbative contributions to the amplitudes and parton distributions. In Sect. IV we show how Basic factorization for 
the quark-hadron amplitudes and quark distributions in hadrons can be reduced to Kt- and Collinear factorizations. 
In Sect. V we focus on a model for the invariant amplitudes T^ u ' > and T ( - s ' 1 . The model is based on description of T^ u ' > 
and T- V in a quasi-resonant way and through the Optical theorem it easily leads to non-perturbative contributions to 
the parton distributions, with expressions of the Breit-Wigner kind both in Basic and in Kt- factorizations. Finally, 
Sect. VI is for concluding remarks. 


II. QUARK-HADRON AMPLITUDES 

In the factorization approach, the quark-hadron amplitudes (QHA) A q are expressed through convolutions of 
perturbative amplitudes A^ pert ^ and non-perturbative amplitudes T as shown in Fig. [3] 

In the Born approximation A^ pert ^ is depicted in Fig. [I] as a one-rung ladder. Adding more ladder rungs to it 
together with inclusion of non-ladder graphs and resumming all such graphs converts the Born amplitude into A^ pert ^. 
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FIG. 3. Factorization of the quark-hadron amplitude. 


In the present paper we do not consider mixing of quark and gluon ladder rungs, i.e. we consider the graphs where 
the vertical quark lines go from the bottom to the top without breaking. 



FIG. 4. Born approximation for the factorized quark-hadron amplitude. 

We begin consideration of the quark-hadron amplitudes A q in Basic factorization, studying the simplest case 
depicted in Fig. U where the perturbative contributions are accounted in the Born approximation and denote such 
distributions B q . In Basic factorization one can use the standard Feynman rules to write down the analytic expression 
corresponding to the graphs in Figs. 3,4. Doing so, we obtain that 


B q = — iAttcXsCf 


d A k u(q)^kT q (k,p)k'j v u(q) 


(2tt) 4 


k 2 k 2 (q + k) 2 




( 6 ) 


where we have used the standard notations: Cf = (iV 2 — l)/(2TV) = 4/3 and a s is the QCD coupling. In Eq. © 
T q corresponds to the lowest blob in Fig. [3] It is altogether non-perturbative object. Throughout the paper we will 
address it as the primary quark-hadron amplitude^. Choosing the Feynman gauge, where d= g^, for the virtual 
gluon and the Sudakov parametrization J3]) for the quark momentum k , we rewrite Eq. m as follows: 


B, = -,^Zw f dad^kA ( “ )l '‘ kT - (k ' V)k ' 1 '‘ U(q) 


87T 3 


k 2 k 2 {q + k) 2 


( 7 ) 


Throughout the paper, for the sake of simplicity, we will treat the external quarks with momentum q as on-shell 
ones, though our reasoning remains valid also when they are off-shell. Introducing the density matrix 


P(p)(q) = ^{<1 + m q )(l - 75 s q ), 


( 8 ) 


1 In Ref. Q non-perturbative contributions to parton distributions in the context of Collinear factorization were called intrinsic contribu¬ 
tions. 
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with q , m q and S q being the quark momentum, mass and spin respectively, we bring Eq. to the following form: 


B„ 


u s Cf 

~ l ~8^r w 


Tr 


dadfid 2 k±- 


p(p)(qhtikT q (k,p)k% 


k 2 k 2 {q + k) 2 


(9) 


We stress that the replacement of Eq. 0 by Eq. @ is not necessary for us but it allows us to carry out a more 
detailed consideration of A q . In particular, we can consider separately the spin-dependent, B q sp "' 1 and independent, 
g(unpoi) q Uar ] £ _i iaiC j ron amplitudes in a simple way: 




dad(3d 2 k±- 


2 (qk)Tr 


kT c 


(unpol ) 


- k 2 Tr 


qTc 


(unpol) 


k 2 k 2 {q + k) 2 


( 10 ) 


B< spm) = ^ 'fm q w J dadpd 2 k± 


2 (S q k) Tr 




(spin) 


- k 2 Tr 


75 S q T q 


(spin) 


k 2 k 2 (q + k ) 2 


( 11 ) 


In Eqs. (1101111) we have replaced the general primary amplitude T q by more specific amplitudes j' q unpol \T q spm \ In 
Eq. (fTTTl) we have neglected a contribution ~ m in p{q) compared to the contribution ~ q. Integrations in Eqs. (11011111 
run over the whole phase space and it is supposed to yield finite results. However, there can be singularities in the 
integrands and they should be regulated. Regulating them with introducing various cut-offs would be unpliysical, 
so the only way out is to impose appropriate constraints on the primary quark-hadron amplitudes j' q unpol \T q spin ' > 
so that to kill the singularities. When the perturbative amplitude A(P ert ) is calculated in the Born approximation, 
the only possible singularities in Eqs. mm are IR singularities at k 2 = 0 and UV singularities which we relate 
to integrations over a. However, when A^ pert ' ) is beyond the Born approximation, there appears another kind of 
singularities called in Ref. fli] rapidity divergences. Below we consider handling these singularities in the framework 
of Basic factorization. 


A. Rapidity divergences of QHA 

Rapidity divergences were investigated first in Ref. i8| and then in Ref. Q in the context of I\t -factorization. 
Detailed investigation of this problem can be found in Ref. 0- In the lowest order of the Perturbative QCD, the 
rapidity divergences come from the graphs in Fig. [5] (and symmetrical graphs as well), where the radiative corrections 
calculated in the first-loop approximation are convoluted with the unintegrated parton distribution <h. Let us stress 
that accumulates both perturbative and non-perturbative corrections. When such convolutions are considered in AY 




FIG. 5. Graphs contributing to rapidity divergences in unintegrated parton distributions. The dashed lines denote cuts. 

-factorization, each of the graphs in Fig. [5] acquires logarithmic divergences arising from integration over momentum 
l + (with l + = (l 0 + l z )/y/ 2). They are called rapidity divergences and they can be got rid of as shown in Ref. @] 
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(when the Feynman gauge is used for the gluon propagators) and then in Ref. {9| for the case of the light-cone gauge. 
In Refs. [$, 9j the rapidity divergences are cured with redefining $. 

Now let us study this situation in Basic factorization. To this end we consider a contribution of the graph in Fig. 6 
to the quark-hadron amplitude in Basic factorization. We remind that there are no cuts in Fig. 6 and the blob T 
accumulates non-perturbative contributions only. One of remarkable features here is that analytic expressions in Basic 



FIG. 6. Graph contributing to quark-hadron amplitude . 

factorization can be obtained by applying standard Feynman rules to the involved graphs. Second important point is 
that one is free to use any gauge for perturbative QCD calculation^ in Basic factorization whereas the blob T in Fig. 6 
is altogether non-perturbative and therefore it is insensitive to the choice of the gauge. Applying the Feynman rules 
to the graph in Fig. [6] and integrating over the loop momentum l, we immediately conclude that this integration yields 
a logarithmic UV-divergent contribution which, being complemented by a similar contribution from the symmetrical 
graph and self-energy graphs, in a conventional way leads to renormalization of the gluon-quark couplings. After 
absorption of such divergent contributions by the couplings, we obtain a renormalized amplitude which is free of 
divergences. Then, applying the Optical theorem to the this construction, we arrive at the parton distributions and 
they are also free of divergences. Obviously, the same treatment can be applied to other UV divergences coming 
from perturbative component A^ pert ^ in higher loops: all of them can be absorbed by renormalizations. Now we focus 
on the divergences resulting from integration of the convolutions in Eqs. (fTUITTD . where the perturbative amplitudes 

Aiper-i) 

are in the Born approximation. 


B. IR and UV stability of QHA 

First of all, let us note that the denominators in Eqs. (110lll|) can become singular in the infra-red (IR) region, 
where k 2 ~ 0. In the case of purely perturbative QCD, IR singularities are conventionally regulated by introducing 
IR cut-offs. In our case there is not any physical reason for that, so we are left with the only way to kill these 
singularities: The primary quark-hadron amplitudes T q should become small at small k 2 : 


j^(unpol) Jji(spin) ^ (/ ^ 2 ) 1+r, (12) 

when k 2 —> 0. Now let us consider the ultra-violet (UV) stability of the convolutions in Eqs. (11011 II ). The integration 
over a in Eqs. (1101111) runs from — oo to oo, so, at large |a| the integrands should decrease fast enough to guarantee 
UV stability. First of all we focus on the integration over a in Eq. m- Taking into consideration that each factor in 
the denominator of Eq. m is ~ a makes that the denominator to be ~ a 3 . The term 2 qk in the numerator depends 
on a because 2 qk = w(/3 — X\a) and the factors k 2 and k are ~ a , which makes 


2 qk k a 2 

k 2 k 2 (q + k) 2 a 3 


(13) 


2 For gauge invariance of Basic factorization see Ref. [jj. 
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This divergence must be regulated by an appropriate decrease of j'( unpo1 ') a t large |a|. The IR stability condition in 
Eq. m states that j'( unpol '> ^ (fc 2 ) 1+I? at small k 2 but it can either disappear or be kept at large |a|. Therefore we 
have two options: 

(A) The factor (fc 2 ) 1+ ’ f survives at large |a|. 

(B) The factor (k 2 ) 1+v disappears at large |a|. 

In the case (A), where IR and UV behaviors of r f^ unpo1 ' 1 are related, r i^' inpol > should behave at large \a\ as follows: 

finnpol) ~ a r,- x = (qjI+ij) [a _1_x ] , (14) 

with x > V > 0. 

IR and UV behaviors of r f( ] unpol > are disconnected in the case (B). It converts Eq. (fUH) into 

f^npoi) ~ a - x _ ( 15 ) 

The first factor in Eq. m corresponds to the term (fc 2 ) 1+,f , while a contribution generating the asymptotic factor in 
the squared brackets has to be specified. We will do it in Sect. V. Now let us consider the spin-dependent amplitudes. 
In order to guarantee their IR stability, the primary spin-dependent amplitude T l q splrl> should also be ~ ( k 2 ) 1+v at 
small k 2 but the situation with its UV stability is more involved than in the unpolarized case. Indeed, the quark spin 
S q can be either in the plane formed by p and q , i.e. S q = Sq, or in the transverse space, where S q = S^~. Depending 

on it, there are the longitudinal spin-dependent amplitude, I3 q and the transverse one, . Now let us consider the 
term 2 mn q S q k in Eq. m for different orientations of the quark spin: When the spin is longitudinal, 

2 m q S q k = 2 m q S^k = w({3 — aqa) (16) 

and k in the trace Tr[kTq\ is also ~ a. In contrast when the spin is transverse, 

2 m q S q k = —2m q {S q k±) (17) 

and therefore S q k does not depend on a. Then, this kj_ should be accompanied by another k± from the trace in 
order to get a non-zero result at integration over the azimuthal angle, i.e. The first term in the numerator of Eq. m 
does not depend on a, while the second term is ~ a. It means that, with T q dropped, the explicit a-dependence of 
A q at large |a| coincides with the one in Eq. (ED: 


and 


Sljk k cf_ 

k 2 k 2 (q + k) 2 a 3 


( 18 ) 


S q k k a 
k 2 k 2 (q + k) 2 a 3 


(19) 


It follows from Eq. (fl8l) states that the a -dependence of the amplitude T q at large |a| is identical to the one of 

rj-\(unpol) 

1 Q 


tW - 

q 

- T( un P oV > „ ( a 1+ri ) [q:- 1 -^] 

(20) 

in the case (A) and 




f'I ~ f^ unpol) ~ (a~ x ) 

(21) 

in the case (B). T f j- can decrease slower: 




T q L - (a 1+? 0 [a" x ] 

(22) 




iii the case (A) and 


T q ~ a 1_x (23) 

in the case (B). Eqs. (1121201 - 1231) guarantee integrability of the convolutions for the quark-hadron amplitudes in Basic 
factorization. These integrability requirements can be used as general theoretical constraints on non-perturbative 
contributions to the amplitudes in Basic factorization (see Ref. Q for detail) and we will use them in the present 
paper. Each of Eqs. (1201221) consists of two factors. The first factor in these equations is universally generated by the 
term (fc 2 ) 1+?) while contributions generating the factors in squared brackets will be specified in Sect. V. 


III. MODELING THE SPINOR STRUCTURE OF T q 

Our next step is to simplify the traces in Eqs. (1101111) . In order to do it, we have to specify the spinor structure of the 
primary QHA T q . By definition, T q is altogether non-perturbative, so specifying its spinor structure can only be done 
on basis of phenomenological considerations. However, any model expression for T q should respect the integrability 
conditions in Eqs. (fl2lll4l2QI22l) . There is the well-known expression for the density matrix of an elementary fermion: 

p(p) = \{P + M ){ 1 - TsS') ~ f ~\(P + M )75 S, (24) 

where M and S are the fermion mass and spin. This expression drives us to approximate T q as follows: 


T q = p T^ u \k 2 , 2 pk) - (p + M) 15 S T^(k 2 ,2pk), (25) 

where p, S are the hadron momentum and spin respectively and T q U \ T q S ^ are scalar functions. Throughout the 
paper we will address them as invariant quark-hadron amplitudes. Substituting T q of Eq. (fl2l) in Eqs. (1101111) and 
calculating the traces, we arrive at the following expressions: 


g(unpol) _ _ 


= —i 


87T 3 

8-7T 3 


dadj3dk\ 


-g C F 


w 


[(9- 


dadpdk\B^ npol \q,k) 


k) 2 - 

k 2 . 




k 2 k 2 


\k 2 k 2 

(rf)(fc 2 ,2pfc)) , 


(rf)(fc 2 ,2pfc)) 


(26) 


were we have denoted Q^ un 'P° 1 ^ the perturbative amplitude in the Born approximation for the forward annihilation of 
unpolarized quark-quark pair. We have neglected contributions ~ Xi t 2 in the numerator of Eq. (1261) and will do it 
in expressions for the spin-dependent amplitudes. These terms, if necessary, can easily be accounted for with more 
accurate implementation of Eq. © to Eqs. (1261) . Let us consider the structure of the integrand in Eq. (1261) in more 
detail. The amplitude in the last brackets is entirely non-perturbative. It is suppose to mimic a transition from 
hadrons to quarks. The fraction in the middle corresponds to the convoluting the perturbative and non-perturbative 
amplitudes. The fraction in the first brackets corresponds to the perturbative amplitude for the forward scattering 
of quarks in the Born approximation. We explicitly wrote the factor le there to remind that this amplitude has the 
s-channel imaginary part. Doing similarly, we obtain an expression for the spin-dependent amplitudes: 



g 2 C F 
1 167r 4 


2 m q Mw 


dad (3d 2 k± 


2{kS q )(kS)-k 2 {S q S) 
k 2 k 2 [(q + k ) 2 + «e] 


rp(S) 

q 


(27) 


Let us consider Eq. m for different orientation of the hadron spin: 

(i) The hadron spin S is in the plane formed by momenta p and q, so for this case we use the notation S = S^. 

(ii) The hadron spin is transverse to this plane. We denote this case as S = S- 1 . 

Amplitude A| for the first case is given by the expression very close to the unpolarized amplitude: 


i?Cll) = 

q 


167T 3 


dad/3dk]^ 


~9 C F 


2mM(SqS^) 
{q + k) 2 


l€ 


\k 2 k 2 


T^\k 2 ,2pk) 


-2^4 J dad/3d 2 k ± B^ ] {q,k) T q l) (k 2 , 2pk), 


(28) 
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whereas the transverse amplitude is given by a different expression: 


B = 


1 


167T 3 

1 


dadftdk\ 


~9 2 C f 


2 mM(S B S B ) 
(q + k) 2 


■ it 




/ dadfddk\(q,k) T q ±] (k 2 ,2pk), 


(29) 


with Bq^\ Bq being the perturbative spin-dependent Born amplitudes. Accounting for perturbative QCD radia¬ 
tive corrections converts the Born amplitudes Bq Unpol \ Bq^\ B^ in Eqs. (1261281291) into perturbative dimensionless 
amplitudes A^ npol \ Aq^, Aq^, remaining the other factors unchanged: 


A^ npol \p,q) = -i± 3 J d p^daA^ ol \q,k) T^ u \k 2 ,2pk), (30) 

A { }\p,q) = -*^3 J d/3-j^daA^\q,k) T^\k 2 ,2pk), 

Ag^&q) = ~ 1 -^3 J dP^^daA^ (q,k) (^f^jT^\k 2 ,2pk). 

Taking the s-imaginary part of Eq. (USD, we arrive at the totally unintegrated, or fully unintegrated as was used in 
Ref. 0, distribution of unpolarized quarks in the hadron £)^ npo1 B in the Born approximation: 

D (un P ° l S) = _i_ J dp ^L da y2 Cf s {p-x- z)\ } (fc 2 ,2 pk) (31) 

= i/ T^ da \p 2c Fd(l-x/p~ z /p)}^^\k 2 ,2pk) 

where x = —q 2 /w, z = k^_/w and is the primary quark distribution of unpolarized quarks in the hadron, 

This object is altogether non-perturbative. Applying the Optical theorem to Eq. (l30l) . we arrive 
at the parton distributions beyond the Born approximation: 


D^ npol \x,q 2 ) 

D^\x,q 2 ) 

D^\x,q 2 ) 


gij J daD(™ p ° l \x/p,q 2 /k 2 )dp('jpj EEf \k 2 ,wa ), 

^2 j ^J^rdaD^(x/(3,q 2 /k 2 ) *M(k 2 ,wa), 

^2 / ^^daD[ B \x/p,q 2 /k 2 )(^^j ^{k 2 ,wa). 


(32) 


IV. REDUCTION OF BASIC FACTORIZATION TO CONVENTIONAL FACTORIZATIONS 

Conventional forms of factorization are Collinear and k F - factorizations. In Ref. Q we described reduction of Basic 
factorization to Kt~ and Collinear factorizations for the Compton scattering amplitudes and DIS structure functions 
without specifying the non-perturbative amplitudes T q . In this Sect, we show that these results perfectly agree 
with our assumption in Eq. (12511 concerning structure of T q . We demonstrate that the parton distributions in both 
conventional factorizations can be obtained with step-by-step reductions of the expressions for ]jk inpor >_ Dq^\ Dq 11 in 
Basic factorization. This reductions are the same for both the parton-hadron amplitudes and parton distributions, 
they are insensitive to spin and stands when the quarks are replaced by gluons. Because of that we consider such 
reductions for a generic parton-hadron distribution D in Basic factorization and skip unessential factors: 


D{x,q 2 ± ) = J d/3^±daD {pert) (x/f3,q 2 ± /k 2 ) ^(wa,k 2 ), 


(33) 
















10 


where D^ pert ^ stands for a perturbative contribution and d/ is the altogether non-perturbative (we address it as 
primary) parton-hadron distribution. Actually, ^(wa, k 2 ) is the starting point for the perturbative evolution. In¬ 
tegration in Eq. (1331) runs over the whole phase space. Let us note that in the literature very often are considered 
purely transverse q: q 2 ~ q\. Because of this reason we will use the notation q\ instead of q 2 in what follows, though 
Eqs. (I32I33D are also valid when q 2 ^ q\. 


A. Reduction to hr -factorization 

In order to reduce Eq. (1331) to kx -factorization, we have to perform integration with respect to a. However, this 
integration should not involve D^ pert \ which, strictly speaking, is impossible because D^ pert ^ depends on k 2 and 
thereby it depends on a\ k 2 = —wa/3 — k\. The only way out is to assume that the main contributions to Eq. (l33l) 
come from the region where 

a < cwx = &i/ (w/3), (34) 

i.e. k 2 k, —k^_. Let us notice that approximating ladder partons virtualities k 2 by their transverse momenta is well- 
known. It is used in all available evolution equations, including DGLAP and BFKL, and now it allows us to convert 
Eq. (1331) into an expression for the unintegrated (transverse momentum dependent[l2|) parton distributions Dkt in 
kx -factorization: 


D K t (x, q±) ~ J ^D { ^p\x/P,q 2 Jk\)^{wP,k\), (35) 

where r = w/3 — q\. Obviously, r ~ w ~ q\ for large x while at very small x one can use r ~ w/3. <f> denotes the 
primary (i.e. non-perturbative) kx -parton distribution. It is related to 4' as follows: 


5>(w/3, k±) 



(36) 


B. Reduction to Collinear factorization 

In Ref. [B] we discussed how to reduce Kx -factorization to Collinear one, using DIS structure functions as an 
example. The same argumentation can be applied to the parton distributions. We briefly repeat it below. In order 
to reduce kx -factorization to Collinear factorization, we should perform integration of Eq. m with respect to k± 
without integrating D^ per " t \ Of course it cannot be done straightforwardly, because D( pert ^ explicitly depends on kj_. 
However, we can do it approximately, assuming a sharp peaked dependence of 't(B/a,Ij) on k\ with maximum at 
k\ = p 1 2 as shown in Fig. [7] The close this dependence is to 8 {k\ — ^t 2 ), the higher is accuracy of the reduction. 
As discussed in Ref. [HJ, the number of such maximums can be unlimited. We remind that d* is non-perturbative, so 
typical values of fi must be of non-perturbative range, n ~ A qcd- After the integration of $ we arrive at Collinear 
factorization convolution: 


D {col) {x,q 2 Jn 2 )^ J ^D < fJ rt) (x//3,q 2 ± / f i 2 )(l)(l3,iJ 2 ), (37) 

with /i being the intrinsic factorization scale and </> being the primary (non-perturbative) integrated parton distribution: 

r 3k 2 

m» 2 )= J Q ^r*(wP,kl), (38) 

where the integration region Q is located around the maximum k ^ = 3 ?■ At the first sight, the form of Collinear 
factorization presented in Eq. (1371) contradicts to the conventional form. Indeed, the scale in Eq. (1371) corresponds 
to the maximum in Fig. 7 and therefore its value is fixed. On the contrary, the conventional form of Collinear 
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FIG. 7. The peaked form of < f>(/3,fc'l) with one maximum 


factorization operates with integrated parton distributions ip(/3,Jl 2 ), where the scale Jl can have any arbitrary value. 
Then, we expect the value of fj, to be of non-perturbative range whereas usually Jl ~ few GeV, i.e. typically Jl^> fi. 
However, this contradiction can easily be solved as was shown in Ref. [g]. The point is that transition from <p(f3,fj, 2 ) 
to ip(/3,Jl 2 ) can be done, applying perturbative evolution in the p 2 -space to 4>(/3,n 2 ) and keeping /3 fixed. It can be 
written symbolically as 


^(/3,fi 2 ) = E(ij, 2 ,ij, 2 ) ® <K/3,Ai 2 ), ( 39 ) 

where E is the evolution operator in the /i 2 -space. Specific expressions for E are different in different perturbative 
approaches (see Appendix E for detail). Eq. rtdflll makes possible to arrive at the conventional unintegrated distribution 
ft 2 ) fixed at an arbitrary scale Jl 2 . In contrast to 4>((3, p 2 ), the distribution tp(/3, ft 2 ) accumulates both perturbative 
and non-perturbative contributions. It is easy to show that our reasoning remains true in the case when $ has several 
maximums or an infinite series of them. This point was discussed in detail in Ref. Q, so we will not do it in the 
present paper. Instead, we focus on modeling invariant amplitudes T^ u,s ^ introduced in Eq. (USD- 


V. MODELING THE INVARIANT QUARK-HADRON AMPLITUDES AND PRIMARY QUARK 

DISTRIBUTIONS 


In this Sect, we suggest a model which mimics non-perturbative QCD contributions in the primary hadron-quark 
invariant amplitudes T^ U ’ S ' ) and in the primary quark distributions in all available forms of factorization. Once again 
we begin with consideration of the invariant amplitudes T^ u,s ^ and then proceed to the quark distributions. 


A. Resonance model for the primary quark-hadron invariant amplitudes 

Amplitudes T^’ 5 ') can be introduced in a model-dependent way only because QCD has not been solved in the 
non-perturbative region. All such models should satisfy several restrictions: 

(i) : The IR stability conditions in Eq. m and the UV stability conditions in Eqs. m-m should be respected because 
they guarantee integrability of the factorization convolutions. We remind that the UV stability conditions derived in 
Sect. II depend on UV-behavior of the factors regulating IR divergences . Namely, Eqs. ()14l20l22|i correspond to the 
case (A) while Eqs. (U5I21I23[) correspond to the case (B). In the present paper we focus on the most UV-divergent 
case (A), although our conclusions hold true for the case (B) as well. 

(ii) : the invariant amplitudes should respect the Optical Theorem, so they should have s-channel imaginary parts. 

(iii) : These amplitudes should guarantee the step-by-step reductions from Basic factorization to kr -factorization and 
then to Collinear factorization described in Sect. IV. 

The expressions for the unpolarized and spin-dependent amplitudes with the longitudinal spin in Eq. m are much 
alike while the expression for the transverse spin amplitude differs from them. Despite this difference, our model 
equally stands for all spin-dependent amplitudes and quark distributions regardless of the spin orientation. To 
describe the invariant primary quark-hadron amplitudes, we suggest a model of the resonance type for Tq U ’ S ' > : 
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Tq U \pk, k 2 ) 
Tq S \pk, k 2 ) 


Ru ( k 2 ) 

{{k - p) 2 - Ml + iT-l) ((k - P y- - Ml + iV 2 ) 
Rs ( k 2 ) 

{{k - pY - Ml + ir 3 ) ((k - pY - Ml + zr 4 ) 


(40) 


where Ru,s ( k 2 ) are supposed to behave as Ru,s ( k 2 ) ~ (k 2 ) 1+v at small k 2 . We need at least two resonances to 
satisfy the UV stability requirement of Eq. m • Indeed, Eq. (l40l) leads to % = 1 while a model with one resonance 
corresponds to y = 0. Formally, Eq. (1401) contains independent parameters Mi t2 , 3 , 4 , and Ti^. 3,4 but we do not see a 
physical reason forbidding to identify Tq J> and Tq S \ which would left us with the parameters M 4) 2 and I \2 only. In 
terms of the Sudakov variables Tq U,S ' > are: 


T^ u \wa,k 2 ) = 
T^Ywot, k 2 ) = 


Ru ( k 2 ) 


(wa + k 2 — p\ + iTi) (wa + k 2 — p 2 + 7 T 2 ) 

_ Rs (fc 2 )_ 

(wa + k 2 — p 3 + *r 3 ) (wa + k 2 — p 2 + *T 4 ) ’ 


(41) 


where k 2 = —wa/3 — k\ and 


p 2 = M 2 -p 2 . (42) 

We suggest that values of p 2 and r 7 - should be within the non-perturbative scale domain, with M 2 > Tj. It is 
convenient to write as the sum of two resonances: 


T^ u \wa,k 2 ) = 


Ru ( k 2 ) 


(p\ — p\) — ?(Ti — r 2 ) [(wa + k 2 — p\ + «r 4 ) (wa + k 2 - p\ + iT 2 )_ 


(43) 


1 1 

(wa + k 2 — p 2 + *r 3 ) (wa + k 2 — p\ + «T 4 ) 

It seems that specifying Ru,s cannot be done unambiguously. We postpone investigating this problem to the future 
while in the present paper we use Ru,s defined as follows: 


T^ s \wa,k 2 ) = 


Rs (M 


(pi - pi) - *( r 3 - r 4 ) 


/ i .2 \ !+v / k 2 \ 1+v 

Rt ' =Al, (fcV^f) ■ R s = ^(pTWJ ’ (44) 

where A u,s and pfj s , (pfj s > 0 ) are independent parameters, though we think that Ru and Rs could coincide which 

would diminish the number of free parameters. It is easy to check now that the expressions for Tq U \ Tq S ^ introduced 
in Eq. (HIT) obey the condition of IR stability in Eq. (fl2l) with arbitrary 77 . In contrast, the value of the UV parameter 
X (introduced in Eqs. (1141201) to guarantee UV stability) is now fixed: x = 1 in Eq. (PHI) . Eqs. (1251411) are suggested for 
invariant amplitudes Tq U ' S ^ in Basic factorization. Reducing Basic factorization to kx -factorization converts 
into new amplitudes Tq U ' S \ They are obtained from Tq U ' S ' 1 by integrating them with respect to a: 


/“®max 

T«(/3,fci)= / daT^(a,k 2 ), 

Jo 

where r = C/, S'. The upper limit of integration, a max should obey Eq. (HU) , so we choose 


(45) 


Omax ~ kj_/(wfJ). 


(46) 
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According to Eq. (l34l) . k 2 ~ 
detail): 

T^ s \p,kl)^\ 

where A Tq U ^ and A Tq S ' ) are 


— k\. The integration leads to the following expression for Tq 1 ' 1 (see Appendix D for 


Ru (fcj.) 

Rs (*i) 


1 


1 


L*i(l - >9)/j9 - A*? + *Tr *1(1 - 0)/P - ^ + *T 2 J 

1 1 


Lfci(l-/ 3 )// 3 -/r| + *r 3 fcl(l-/ 9 )/ 0 -M 4 + *T 4 J 


■A 


+ ATf, 


(47) 


1 


fcj + /Ltf — iTi 


AT (c/) = ---In , _ , 

9 {41 - Mi) + *(Tr - r 2 ) V fjiIpiI ’ 

A f(S) = _I_ hl f kl+ti-tTa ' 

9 (/il-^) + *(r 3 -r 4 ) + 


(48) 


They depend on k± very slowly and they can be neglected at large k\. 


B. Primary quark distributions 


The Optical theorem relates the s-channel imaginary parts of T^ ,s ) an d r f(u,S) ^g p r i mar y quark distributions 
^u,S in Basic factorization and to unintegrated (or ) quark distributions <&u,S in fcr -factorization respectively. So 
applying the Optical theorem, we obtain the following expression for the primary quark distribution 'k r in Basic 
factorization: 


T;/(to, k 2 ) = 


_ 1 Ru (k 2 ) 


n {41 - 4b 


Ti 


4's(u'a, k 2 ) = — 


1 Rs{k 2 ) 


n (41 - 41) 


(■wa + k 2 — 41 ) + Tf (to + k 2 — 4 b + T 2 

_Tg_ U _ 

(to + k 2 - 4lf + r§ (to + k 2 - 41? + r 2 


and a similar expression for the primary quark distribution $>u,S in kx -factorization: 


(49) 


$u(P, k 2 ± ) = —Ru (& 1 ) 


$s(P,ki) = -R s (fci) 


+ 


(*1( i - P)/P - 4i? + r? (*i(i - pyp - 4i? + r| 

r 4 


(*i(i - P)/P - 41? + + (fci(i - /3)//? - 41Y + r| 


(50) 


Obviously, the expressions in Eqs. (1491501) are of the Breit-Wigner type. Substituting Eq. (l50l) in Eq. (1351) and 
integrating over k^_, we arrive at the quark parton distribution Dj Co1 ^ in kx -factorization, where the non-perturbative 
contributions i.e. the unintegrated parton distributions are specified: 


D^\x, q D = i f J f d -^&r ,) (x/Dyjdi)R U (ty 


(51) 


Rs (k 2 ± ) 


(fci(i - p)//3 - 41) + r? {kl (i pyp - 4b + r| 

Df T \x,qi) = -\ X % f ^D^ rt \x/P,k 2 j q 2 ± 

71 JX P JO K ± 

r 3 , r 4 


+ 


(*I( i - P)/P - 41) + r§ (fci(i - pyp - 41) + r 4 
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where r is defined in Eq. (l35l) . Let us consider the k± -dependence in Eqs. (1501351) in more detail. Obviously, the 
structures of expressions for D^ T ' 1 and D < g T ' > (or and 4 ) 5 ) are quite similar, so we consider D^ T ' 1 only. Then, the 
expression in the squared brackets in Eq. fl35| , i.e. $(7 of Eq. (1501) . is symmetric with respect to replacement 1^2. 
Each term in the parentheses has a peaked form, with maximums at k\ = /if 2 . The less IT^, the sharper the peaks 
are. We remind that Ru,s ~ {k^_) 1+v at small k\. By definition, see Eq. (H21) . /if 2 = ^12 — p 2 , so they can be either 

positive or negative while k\ cannot be negative. In any case the both terms in $(7 and $s contribute to D^ T g but a 
result of interference of the two peaks depends on values of the parameters. There are possible three particular cases: 
Case (i): both /if and /if are positive. 

In this the both maximums are within the integration region of Eq. (1351) and interference of the two peaks generates 
various forms of $[/(/3, k\) ranging from the picture with two isolated peaks to a kind of plateau, depending on values 

of r lj2 . 

Case (ii): /if > 0 and /if < 0 or vice versa. 

Here the peak from the first term in Eq. (l50l) combines with a tail of the contribution of the second term whose 
maximum is beyond the integration region of Eq. (1351) . The resulting picture has a resemblance to the dual model 
combining a resonant and a constant term. 

Case (iii): both /if and /if are negative. 

The both maximums now are out of the integration region, so tails of the peaks, taken by themselves, generate a 
form slow decreasing with growth of k\. However, this slope is affected by an impact of Rjj- We remind that Ru = 0 
at = 0 . 


C. Primary quark distributions in Collinear factorization 

Performing integration over k\ in Eq. (1351) . we arrive at the parton distributions D^ ol ' > in Collinear factorization. 
Presuming that parameter is small, we write the result of the integration in the following form (see Appendix C 
for detail): 


D { u° l) (x, q i) 


^-D^ ert) (x/P, q 2 j_/nl)<l>u(P, Mi) + J ^-D^ ert) (x/p, q 2 J fJ%)<t>u{P, m!), 


(52) 


with 


4>u(P,P 1 ) 

2 ) 


1 f dk\ Ru (fci) Ti 

* Jth (fci(l-/?)//3-M?) 2 +lf 
1 f dk\ Ru (fcj) r 2 

*Jn 2 k 2 ± {kl{l-P)/p-tif+Tl' 


(53) 


where the integration regions STi = 0 ) f)[ 0 , w] and STi = il' 2 f][ 0 , w], with the subregions Of, Of being located around 
the maximums of the peaks. Formally, the both terms in Eq. (1531) contribute to 4>u at any signs of /if,/if, but in the 
limit of sharp peaks these contributions have different weights: At /if > 0, /if > 0 the both terms contribute equally: 


tu » Ru {p\P/{ 1 - P)) /Ml + Ru {iAP/ (1 - P) /m 2 + o(r 1 ,r 2 ). (54) 

Mostly the first term contributes, when /if > 0, /if < 0: 

Ru(p\P/{ 1-P)) /pl + OiTj. (55) 

and vice versa. Finally, at /if, /if < 0 only tails of the both peaks contribute and therefore 4>u is small and flat 
compared to the previous cases: 


cf>u ~ const. 


(56) 


(see Appendix E for detail). When /if > 0,/if > 0 
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Vuiu, M 2 ) = J j- P u [E(w, /x 2 , ii\)Ru ( M f/3/(l - P)) M 2 ” 2 + E(w, n 2 , &)Ru (mI/W - P)) M 2 “ 2 ] , (57) 

Ru (n*p/( 1 - /?)) nf = Xup. (58) 

Combining Eqs. <E3> and USED, integrating over /? and remembering that at small x essential values of ui are small 
leads to the following expression for (£[/(w,/x 2 ) (see Appendix E for detail): 


<Pu{u,H 2 ) = J y/ 3“ +1 


^E(u, h 2 , h\) + ^-E(u, [i 2 , nl) 
Mi Mi 


—jE(lv, /x 2 , /x?) + ^E(w, /x 2 , M 2 )- 

Mi M 2 


(59) 


VI. CONCLUSION 

In the present paper we have considered the quark-hadron scattering amplitudes and distributions of polarized 
and unpolarized quarks in hadrons in the framework of the factorization concept where the both amplitudes and 
distributions are expressed through convolutions of the perturbative and non-perturbative components. We began 
with considering the quark-hadron amplitudes in Basic factorization where integration over momenta of connecting 
partons runs over the whole phase space and obtained the conditions for the factorization convolution to be stable both 
in IR and UV regions. Then we demonstrated how to reduce Basic factorization to Kt- and Collinear factorizations. 
We suggested a Resonance Model for non-perturbative contributions to the unpolarized and spin-dependent parton- 
hadron scattering amplitudes. This model is based on the simple argumentation: after emitting an active quark by a 
hadron, the remaining colored quark-gluon state cannot be stable and therefore it can be described by quasi-resonant 
expressions. We needed at least two resonances in Basic factorization and this remained true when Basic factorization 
was reduced to Kt -factorization. Applying the Optical theorem to the Resonance Model provided us first with the 
expressions of the Breit-Wigner type for non-perturbative (primary) contributions to the quark distributions in Basic 
and Kt -factorizations and then, after one more reduction, to the parton distributions in Collinear factorization. To 
conclude, let us notice that the Resonance Model can also be used for analysis of the non-singlet components of the 
DIS structure functions. 
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Appendix A: Amplitude for the forward Compton scattering off a gluon in the box approximation 


e 2 a /' 

M^xp = / d/3'da'dk' 2 


yi/fW _i_ m ( 2 ) _i_ 

pvXp ' pbvXp ' fiuXp 


KX = 


KX = 


Tr 


(q + k ,S j 'Yph'^yx {k' - kj jpk' 


Tr 


k' 2 k' 2 (q + k') 2 {k' -k) 2 
(q + k'^j ^pk'^p (k’ + k"j 7 \k' 


KX = 


Tr 


k' 2 k' 2 (q + k') 2 (k' + k) 2 
lv (q + k' - fc) 7 P (k' - k \ hik'lx [k' + qj 


k' 2 {q + k') 2 {k' — k) 2 (q + k' — k ) 2 


(Al) 


(A2) 
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Appendix B: Projection operators for forward Compton amplitudes 

The conventional of dealing with the forward Compton scattering amplitude A is, in the first place, to simplify 
their tensor structure. To this end, Ap V is represented as an expansion of A pv into the series of simpler tensors, each 
multiplied by an invariant amplitude. Such tensors are called projection operators. Through the Optical theorem the 
invariant amplitudes are related to the DIS structure functions. 

In the case of the unpolarized Compton scattering such an expansion looks as follows: 

V = P ^ A 1 + P ff A 2, (Bl) 


where 


p< in) = ~g^ + p$ = (i/w) (pp. - q^{pq/q 2 )) {p v - q»(pq/q 2 )) (B2) 

are the projection operators and A\ 1 A 2 are invariant amplitudes. According to the Optical theorem 

Fi = -9h4i, F 2 = -9fA 2 . (B3) 

7 r 7r 

Similarly, for the polarized Compton scattering 

V = P $ A 3 + P $ A ^ ( B4 ) 

where 

P$ = le^xpMqxSp, P A) = te^xpMqx [S p - p p {qS/qp)] , (B5) 

with M and S being the hadron mass and spin respectively, and A 3i4 are spin-dependent invariant amplitudes. The 
Optical theorem states that 


9i = -9fA 3 , 32 = -SA 4 . (B6) 

7 r 7 r 

All operators Pffl respect the electromagnetic current conservation: q t ,P\lu = q„Pffl = 0. It is convenient to 
introduce the longitudinal, 5^ and transverse, S' _L components of the spin, so that S ± p = S ± q = 0 and = p p (qS/pq). 
In such terms Eq. (1B41) can be written as follows: 


Apu — ic pu\pAl qx 


SpA 3 + (A 3 + A 4 ) 


= itpvXpMqx 


+ S^A 1 - 


(B7) 


This expression is useful for practical attributing different terms in the spin-dependent A pu to proper invariant 
amplitudes. In the unpolarized case one can use the simple rule: expressions ~ g pl , contribute to Ai while expressions 
~ PpPu/pq form A 2 . In contrast, the gauge invariance admits adding arbitrary terms ~ q p , q v . 


Appendix C: Convolutions involving the Breit-Wigner formula 

Let us consider the following convolution: 



7T 


/ oo 

dxf{x)-~ 

-00 (A 


r 

z 0 ) 2 + r 2 ‘ 


(Cl) 


Replacing x by t, with t = (x — Xq )/T, we convert Eq. (IC1I) into 
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i r°° r 

F = - dtf(tT + Xq)- 

7T I 


t 2 + 1' 

At small r, we can expand f(tT + xo) in the power series and retain several terms: 

f(tT + z 0 ) = f(x 0 ) + f'(xo)tT + o(r 2 ). 
Substituting Eq. (Oil in Oil and integrating m yields 


(C2) 


(C3) 


F = f(x 0 ) + O(T). 

The first term in Eq. (IC4[) corresponds to the well-known representation of the 5 -function: 


(C4) 


- lim 9 € = 6(x) 

7 r £->o x 1 + 


(C5) 


Appendix D: Integration in Eq. (1451) 

A generic expression to integrate can be written as 


T = 



{a 


da 
A) (a 


B) 


1 

W~B) 





(A-B) [ 


In 


- A 






(Dl) 

Assuming that a max A, B \ A — B\ allows us to expand the logarithms into the power series and retain the first 
terms only: 


T Ki- 
2 


A a max B 


1 


(A-B) 


In (A/B). 


(D2) 


We have written Eq. m in the symmetrical form with respect to A, B because Eq. m has this feature. 


Appendix E: Evolving the factorization scale in Collinear factorization 


Using the Mellin transform, we can rewrite Eq. m as follows: 

/ lo ° duj 

—x~ u Cu(uj)E(uj,q 2 : iJ, 2 )(l)(uj,iJ, 2 ), (El) 

-zoo 27TZ 

where the primary quark distribution is fixed at the scale /x, with fj, 2 < q\. E(uj, q 2 ^, fj 2 ) is a generic notation for an 
operator evolving the distribution </>(w,/z 2 ) from the factorization scale [F to q\, while C(u>) is responsible for the x 
-evolution. Choosing a scale /x such that fi 2 <J1 < q 2 ^ and representing E(oj, q\, /z 2 ) as 


F(i^,q 2 j_,E 2 ) = E(u},q 2 ± ^i 2 )E(uj,^ 2 ,fj, 2 ), 


(E2) 


we bring £)( co b to the conventional form 



dui 

2m 


Cjj(u})E(oj, q 2 ,iJ, 2 )(fi(uj, n 2 ), 


(E3) 


where 


^(w,At 2 ) = E(u,n 2 ,ii 2 )(t)(l3,n 2 ), 


(E4) 
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which corresponds to Eq. (1391) . Actually, tp(aj,Jl 2 ) is the conventional parton distribution in the ui -space (momentum 
space). It is fixed at an arbitrary scale ju 2 and related to the standard integrated distribution Sq(x,Jl 2 ) by the Mellin 
transform. The evolution operator E(ui, q 2 , /r 2 ) is expressed in different terms, depending on the perturbative approach 
in use. For instance, in LO DGLAP with fixed a s it is given by 


E = exp [a s 7 o(w) In (q 2 /fJ?)\ , (E5) 

with 70 being the LO DGLAP anomalous dimension and Cf = 4/3. When in LO DGLAP a s is running and the 
standard parametrization a s = a s use< 4 i n the Feynman graphs, Eq. (IE5I) is changed by 


E = 


jVy ' o/6 


(E6) 


with b being the first coefficient of the /3 -function. The parametrization a s = a s (fcj_) should not be used at small x 
(see Ref. [14j). When it is replied by the appropriate parametrization and when the total resummation of the leading 
logarithms is done, Eq. (IE 6 I) is replaced by 


E = exp [h(uj) ln(g 2 //x 2 )] , (E7) 

where h(oj) is a new anomalous dimension. It accounts for the total resummation of the leading double-logarithmic 
contributions and running QCD coupling effects (see Ref. [ll| and overview 0 for detail). 
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